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Abstract 

Quantum chromodynamics (QCD) contains field configurations which can be char- 
acterized by a topological invariant, the winding number Q^. Configurations with 
nonzero break the charge-parity {CV) symmetry of QCD. We consider a novel 
mechanism by which these configurations can separate charge in the presence of 
a background magnetic field - the "Chiral Magnetic Effect". We argue that suffi- 
ciently large magnetic fields are created in heavy ion collisions so that the Chiral 
Magnetic Effect causes preferential emission of charged particles along the direction 
of angular momentum. Since separation of charge is CT^-odd, any observation of 
the Chiral Magnetic Effect could provide a clear demonstration of the topological 
nature of the QCD vacuum. We give an estimate of the effect and conclude that it 
might be observed experimentally. 



1 Introduction 



A remarkable property of gauge theories is the existence of topologically non- 
trivial configurations of gauge fields [Ij. In Minkowski space-time, these con- 
figurations are associated with tunneling between different states which are 
related by topologically non-trivial gauge transformations [2lf3j. These differ- 
ent states are characterized by the Chern-Simons topological charge |3]. Tran- 
sitions between these different states which involve changes in the topological 
charge induce anomalous processes, that is, processes which are forbidden by 
the classical equations of motion of the theory, but when the theory is quan- 
tized, become allowed [5|6] . The amplitudes of such transitions vanish order 
by order in perturbation theory, and are associated with intrinsically non- 
perturbative phenomena involving many quanta. In electroweak theory, such 
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transitions are associated with B + L violation, where B is baryon number and 
L is lepton number. In QCD, they are tied to chirahty non-conservation. Chi- 
rahty conservation is an exact symmetry of the classical equations of motion 
of QCD and of perturbation theory in the zero quark mass limit. The zero 
quark mass limit should be a good approximation for high energy processes 

It was originally thought that non-perturbative phenomena associated with 
the underlying topology were suppressed by factors such as e~^'^/" where a is 
the interaction strength of the gauge theory. This is typical of results from 
computation of quantum tunneling between states with different topological 
charge. It was a surprise when Manton and Klinkhamer proposed that for 
processes at high temperature, that the exponential suppression might disap- 
pear [7l8|. This is because at high temperatures, there is sufficient energy to 
pass classically over the barrier which separates states of different topological 
charge. Kuzmin, Rubakov and Shaposhnikov argued that the rate would be 
sufficiently large so that these processes might be important for electroweak 
baryogenesis PITO] . The rates for such processes were computed and it 
was eventually understood that such rates were not exponentially suppressed 
because they involved the coherent interaction of many low-energy quanta in 
an initial state going to many quanta in a final state ^i2j . Such rates should be 
large therefore not only at finite temperature but also in out-of-equilibrium 
processes which involve many particles. It was also shown that the computa- 
tions for electroweak baryon number violation have a counterpart in helicity 
non-conservation in QCD [T3] . 

An interesting feature of the transitions changing topological charge is that 
they involve V- and CP-odd field configurations. Because of this, V- and CV- 
invariance of strong interactions is no longer "natural". These symmetries 
can however be imposed dynamically by postulating the existence of axions 
[T^15|16] : see [T7] for reviews. In any case, V- and CP-invariance of strong 
interactions has to be related to the vacuum structure. It is thus possible that 
an excited vacuum domain which may be produced in heavy ion collisions as 
proposed by Lee and Wick [18] can break parity and and CV spontaneously 
[T9] . Some consequences of such a scenario had been studied by Morley and 
Schmidt ^20j, who proposed as a signature CP-odd correlations between the 
spins of the protons produced in a heavy ion collision. Other realizations of 
an excited vacuum domain (including the ones with parity-odd fiuctuations) 
have been considered in Refs |21j . 

It was proposed that in the vicinity of the deconfinement phase transition 
QCD vacuum can possess metastable domains leading to P and CP violation 
[22] ■ It was also suggested that this phenomenon would manifest itself in spe- 
cific correlations of pion momenta [22|23] . (For related studies of metastable 
vacuum states, especially in supersymmetric theories, see [2^25f26j ). The 
existence of P*-odd bubbles had been inferred from the analysis of an effec- 
tive chiral theory incorporating axial anomaly [22] ; such bubbles may also be 
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viewed as the space-time domains occupied by gauge field configurations with 
non-trivial winding number. 

The existence of metastable V-odd bubbles does not contradict the Vafa- 
Witten theorem [27] stating that V and CV cannot be broken in the true 
ground state of QCD for = 0. Moreover, this theorem does not apply to 
QCD matter at finite isospin density [28] and finite temperature [29], where 
Lorentz-non-invariant V-odd operators are allowed to have non-zero expecta- 
tion values. For matter at zero temperature and finite baryon density, Migdal 
[50] long time ago advocated the possibility of parity-breaking pion conden- 
sate phase (for a recent work on P-violation in cold dense matter, see [31] 
and references therein). Degenerate vacuum states with opposite parity were 
found [32] in the superconducting phase of QCD. A parity broken phase also 
exists in lattice QCD with Wilson fermions [33], but this phenomenon has 
been recognized as a lattice artifact for the case of mass-degenerate quarks; 
spontaneous V and CV breaking similar to the Dashen's phenomenon p3] 
can however occur for non-physical values of quark masses [35]. P-even, but 
C-odd metastable states have also been argued to exist in hot gauge theories 
[35] . The conditions for the applicability of Vafa-Witten theorem have been 
repeatedly re-examined in recent years [37] • 

Several dynamical scenarios for the decay of V-odd bubbles have been 
considered [38] , and numerical lattice calculations of the fluctuations of topo- 
logical charge in classical Yang-Mills fields have been performed [39lH0] . A 
closely related topic is the role of topological effects in high-energy scattering 
and multi-particle production which has been studied in Refs [^42ll43f44f45j . 
The studies of V- and CP-odd correlations of pion momenta [^6|7F] have 
shown that such measurements are in principle feasible but would require 
large event samples. 

Some time ago, it was pointed out that the presence of non-zero angular 
momentum (or equivalently of magnetic field) in heavy ion collisions makes 
it possible for P- and CP-odd domains to induce charge separation in the 
produced particles [IB] . This is because the P- and CP-breaking term carries 
net chirality and thus generates asymmetry between left- and right-handed 
fermions ~ "ipLipR — "^RipL 7^ 0. In the presence of a magnetic field a P- 
and CP-odd domain can generate chirality not by flipping the spins of the 
quarks, but by inducing up-down asymmetry (as measured with respect to 
the symmetry axis of the angular momentum) in the production of quarks 
and anti-quarks [H]. The magnitude of the effect has been estimated, and 
found to be within experimental reach. A very useful experimental observ- 
able measuring the charge separation with respect to reaction plane (which 
is perpendicular to the angular momentum and magnetic field axis) has been 
proposed by Voloshin [IH] shortly afterwards. The first result of a dedicated 
experimental study of the effect has been recently presented by the STAR 
Collaboration at RHIC [50]. 

In Ref. [51] it has been noticed that the effect of charge separation can also 
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be understood in the following way. At finite ^-angle (as present effectively in a 
V- and CP-odd vacuum domain), a magnetic field B induces through the axial 
anomaly a parallel electric field ^ 6 E, and a corresponding electric dipole 
moment. One may consider as a theoretical illustration of this phenomenon 
the transformation, at finite 6, of a magnetic monopole into a "dyon" with 
an electric charge ~ — e 6/27i demonstrated by Witten [52]. In an off-central 
heavy ion collision, a V- and CP-odd vacuum domain was shown [5T] to 
induce an electric field perpendicular to the reaction plane of magnitude ~ 
— (e 9/2tt) /, where / ~ 6 is the angular momentum in a collision at impact 
parameter h. The corresponding electric charge separation has been estimated 
and found amenable to experimental observation. 

In this paper we discuss a novel mechanism for charge separation. The 
topological charge changing transitions provide the P- and CP- violation 
necessary for charge separation. The variance of the net topological charge 
change is proportional to the total number of topological charge changing 
transitions. Hence if sufficiently hot matter is created in heavy ion collisions 
so that topological charge transitions can take place, we expect on average in 
each event a finite amount of topological charge change. 

Charge separation needs a symmetry axis along which the separation can 
take place. The only symmetry axis in a heavy ion collision is angular mo- 
mentum which points in the direction perpendicular to the reaction plane. In 
central collisions there is no symmetry axis, so in that case charge separation 
should vanish. 

If charge separation would manifest itself mainly into pions it should vio- 
late isospin symmetry (interchange of up and down quarks). This is because 
under an isospin transformation a vr"*" becomes a vr" meson, so in an isospin 
symmetric situation vr"*" should behave similar to vr^, hence in that case they 
can never separate. Isospin violation can occur because up and down quarks 
have a different mass, because our initial state could be isospin asymmetric, 
and because up and down quarks have different charge, or because of the 
presence of an electromagnetic field. 

Since QCD contains configurations that break the CP symmetry, but not 
the C symmetry, the charge separation mechanism we are looking for should 
not violate C. However, a particle with positive charge transforms to a par- 
ticle with negative charge under charge conjugation. In order to maintain C 
invariance, some important ingredient of the mechanism that separates charge 
should transform under C as well. The only possibility lies in the initial state, 
since the incoming ions have a positive electric charge, they transform under 
C. So our mechanism should somehow depend on the electric charges of the 
ions. 

Summarizing, the charge separation mechanism should provide a symme- 
try axis, it should violate isospin and depend on the charges of the colliding 
ions. Therefore it is very natural that the (electromagnetic) magnetic field 
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is crucial for charge separation. This magnetic field depends on the electric 
charges of the colliding nuclei, so that C invariance is maintained. The mag- 
netic field provides a symmetry axis; midway between the two colliding nuclei, 
the magnetic field is parallel to the angular momentum vector. Moreover the 
magnetic field causes violation of isospin since up and down quarks will be 
polarized oppositely in a magnetic field. 

We can understand the effect we will discuss in this article in a simple 
example. Suppose wc have a box with some net topological charge. This 
means that due to the anomaly, there is a net helicity for the quarks inside the 
box. If the box is at high enough temperature, then perturbative interactions 
which involve spin flip are negligible, since T S> mquark, so the net helicity is 
directly proportional to the net topological charge. Now imagine we apply an 
cicctrodynamic magnetic field. This will orient the spin of positively charged 
quarks (and anitquarks) along the direction of the magnetic field, due to the 
magnetic moment interaction. Likewise, the spins of negatively charged quarks 
(and anitquarks) will become oriented opposite to the direction of magnetic 
field. Now because there is net helicity, the orientation of spin is preferentially 
in the direction of momentum of the particles. Positively charged quarks and 
antiquarks move in the direction of the magnetic field, and negatively charged 
quarks and antiquarks move opposite to it. Therefore an electromagnetic 
current is set up in the direction of the magnetic field. We therefore call 
this the Chiral Magnetic Effect. The Chiral Magnetic Effect induces a vector 
current. For systems of a finite size, this current results in charge separation 
or an induced electric dipole moment. 

The separation is as required C even, since under C both the magnetic field 
and the vector current change sign. Under P, the current changes sign, but 
the magnetic field does not, so it is V odd, and CV odd. 

This article is organized as follows. In Section 2 and 3 we review the axial 
anomaly and the QCD vacuum structure respectively. In Section 4 we discuss 
in detail the Chiral Magnetic Effect. Section 5 is devoted to the consequences 
of this effect in heavy ion collisions. In Section 6 we define some observables 
and Section 7 contains a discussion of our results. The derivation of the 
magnetic field is discussed in Appendix A. An alternative derivation of the 
Chiral Magnetic Effect is given in Appendix B. 



2 The Axial Anomaly 

In this section, we review some properties of topologically non-trivial gauge 
fields. We also review the relationship between topology and anomalous pro- 
cesses. 

All gauge field configurations which have finite action can be characterized 
by a topological invariant, the winding number Q^. The winding number is 



5 



an integer and given by 

Qw = d'xF^^^PreZ. (1) 

Here g denotes the QCD couphng constant with generators normahzed as 
tr tah = Sabf^- The gluonic field tensor and its dual are respectively F^^, and 
~ h^ij.u'"^ ^pa- Configurations with nonzero Qw lead to non-conservation 
of axial currents even in the chiral limit. This can be seen from the axial 
Ward-identity which reads 

d^jl = 2 mSf^l^^f)A - ^C^r, (2) 

/ 

here Nf denotes the number of quark flavors, V"/ quark field, and mf the 
mass of a quark. The axial current is given by 

il = Y{i'n,^l5i^f)A- (3) 
/ 

We write {0)a to denote an average over fermionic fields only in the back- 
ground of a particular gauge field configuration A^. To obtain the full expec- 
tation value (O) one still has to perform the integration over the gluonic fields 
A^. In this article we will only consider the chiral limit, i.e. = 0. 

Assuming that initially at t = — oo we have an equal number of right- 
handed and left-handed fermion|3, i-e. Nji = N^, it follows from the axial 
Ward identity that at t = oo 

{Nl - N^)t=oo = ^NfQ^. (4) 

This shows that configurations with positive Qw convert right-handed fermions 
into left-handed ones. 

Let us now write the number of right- or left-handed fermions of a par- 
ticular fiavor as N^j^. So Nji^^ = YIj-^rl- Assuming that initially for each 
fiavor separately = N[, one can show that t = oo 

{Nl - N]^)t=oo = {Nt - N'^)t=oo (5) 

where in this equation u and d can denote any quark fiavor. This implies that 
in the chiral limit, the number of converted up and down quarks is equal. 

^ With fermions we mean quarks and antiquarks. We write -/V^ for the total number 
of right-handed fermions, which is the sum of the number of right-handed particles and 
right-handed antiparticles. Here right-handed massless fermions have spin and momentum 
parallel, while left-handed massless fermions have spin and momentum antiparallel. The 
volume integral over V'Jj'^fl however, is equal to the total number of right-handed particles 
minus the total number of left-handed antiparticles. 
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3 The QCD Vacuum 



In the classical vacuum of QCD the gauge field has to be a pure gauge so that 
the energy density is minimal. In the temporal gauge {Aq = 0) this implies 
Ai{x) = ^U{x)diW (x), where U{x) is an element of the gauge group SU(3). 
The different classical vacua can be characterized by a topological invariant, 
the winding number which is an integer and given by 

n„ = d^xe^^'^tr [iU^d,U)iU^d,U)iU^dkU)] . (6) 

Now one can show that if a gauge field configuration goes to a pure gauge at 
infinity and has nonzero Qw it induces a transition from one classical vacuum 
to another, more precisely 

Qw = ny,{t = oo) - n„{t = -oo). (7) 

At zero temperature such transition requires tunneling through a potential 
barrier, which will suppress the transition rate exponentially. In this case 
the main contribution to the transition rate comes from fiuctuations around 
instantons (which are minima of the classical Euclidean action). It was found 
by 't Hooft [2] (see also [53]) in QCD without fermions that the number of 
transitions A^^'' with Q„ = ±1 per unit volume and time equals 



,,,,, f2nY / 27r\ 1 
* 0.0015 — exp ' ' 



5' 



d^xdtdp \Q^s7 V OLs J P' 

here p denotes the size of an instanton and as = g'^{p)/ (47r) the renormalized 
coupling constant. 

At finite temperature the size of an instanton cannot be larger than 1 / {gT) 
because of Debye screening. Because large instantons cannot contribute to the 
transition rate at large T anymore and g{T) goes to zero due to asymptotic 
freedom, the transition rate due to instantons will go down if the temperature 
is raised. More specifically Pisarski and Yaffe [51] computed the transition 
rate in the absence of fermions and found 



"^^^ - exp(-27rVT2 - 18A(7rpT)) -^4V ' (9) 



d?x dt dp d?x dt dp 



T=0 



where A{X) is a function which behaves as — log(A)/6 for A — >• cxd and as 
-AV36 for A [55]. 

Hence for high temperatures the transition rate due to instantons becomes 
extremely small. This also happens in the electroweak theory where the cou- 
pling constant is always tiny. In this case transitions between different vacua 
lead to violation of baryon plus lepton {B + L) number [2]. However, it was 
realized that the electroweak theory contains static solutions which have finite 
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energy and half integer winding number [8] . Such a solution with minimal en- 
ergy is called a sphaleron (ready to fall). If the temperature is higher than 
the energy of the sphaleron it is likely that one can go over the barrier instead 
of tunneling. This increases the transition rate and hence the rate for baryon 
number violation enormously Pll0|ll|12j . 

It was realized that these sphaleron configurations also exist in QCD [13J . 
The typical energy of the sphalerons in QCD is Aqcd ~ 200 MeV. Hence also 
in QCD at high temperatures one can go over the barrier which leads to an 
enormous increase in the transition rate. According to Refs. ^56j-^57j the rate 
for SU(2) Yang-Mills theory is 

"^^^ ~25.4a^T^ (10) 



d^a; dt 



The rate has not been directly measured for QCD, but if we assume A''^ scaling, 
we estimate the rate for QCD to be 

4^ = T^^ 192.8a|T^ (11) 

There is a subtlety in these considerations for massless quarks. For massless 
quarks, the overlap between a state with zero Chern-Simons charge, and one 
with one unit of charge vanishes, 

<n^ = 1 \ = >= 0, (12) 

since in the path integral representation for such a transition, the fermion 
determinant has a zero mode and vanishes. However, there is a non-zero 
contribution for 

< = 1 I ipip \ = > (13) 

since in the path integral representation for this quantity, the zero from the 
determinant is canceled by a singularity from the quark propagator which 
represents the operator ijjip. The overall expression is proportional to zero 
mode wavefunctions. This operator corresponds to a chirality flip, and as it 
should satisfies the anomaly relation between the change of helicity and the 
change of Chern-Simons charge. 

The total rate of transition is of course the sum of the rates of the lowering 
and raising transitions, 

^ = V (14) 

d^xdt ^d^xdt ^ ' 



4 The Chiral Magnetic Effect 

We will show in this section that gauge field configurations with nonzero Qw 
can separate charge in the presence of a background (electromagnetic) mag- 
netic field B. This we call the Chiral Magnetic Effect. We consider the 
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most ideal situation in which all quark masses vanish. This should be a good 
approximation in the deconfined, chiral symmetry restored phase. First we 
consider the case with a very large magnetic field, which means for us that 
eB is much larger than the momentum squared of all particles, i.e. eB ^ p^. 
Then we will study the effect of a moderate field. 

4-1 Large magnetic field 

Initially we assume that we have a very small number of left and right-handed 
fermions. This is because our initial state has little fluctuation in the helic- 
ities of particles. This is a reasonable assumption for an initial state typical 
of a Color Glass Condensate, where the coupling is weak and the fluctua- 
tions should be generated by quantum mechanical tunneling, and therefore 
suppressed as exp(— 27r/a5) [55]|59f42] . At later time, sphaleron-like transi- 
tions should be possible since the system becomes an ensemble of classical 
conflgurations with energy typically greater than the barrier [60] . 

Since the magnetic fleld is so large, all particles will be found in the lowest 
Landau level. Hence their spin is aligned along the magnetic fleld and they 
can also only move along the magnetic fleld. Quarks with opposite charges 
have their spins aligned in different directions. If B points in the z-direction 
positively charged right-handed fermions and negatively charged left-handed 
fermions will move in the positive z-direction (upwards). At the same time 
positively charged left-handed fermions and negatively charged right-handed 
fermions will move downwards. We illustrated this situation in Fig. [H 

The fermions will then interact with a gauge fleld conflguration, so that 
some of them will change their helicity. After the interaction with the gauge 
fleld conflguration has taken place we flnd that (since we assumed the chiral 
limit) 

{N[ - N^) = 2Q^, (15) 

so that if is nonzero we obtain a difference between the number of left 
and right-handed fermions that is the same for each flavor. The fermions can 
only change helicity by reversing their momenta, since spin flip is energeti- 
cally suppressed in a large magnetic fleld (we assumed eB ^ p^). So after 
the interaction has taken place all positively charged right-handed fermions 
and negatively charged left-handed fermions will still move upwards. At the 
same time positively charged left-handed quarks and negatively charged right- 
handed quarks will still move downwards. However, now there is a difference 
between the number of right-handed and left-handed fermions. As a result an 
electromagnetic current is generated along the direction oi B. If we assume 
that initially all fermions are within a flnite volume, this current will induce 
a charge difference Qe between opposite sides a plane perpendicular to the 
magnetic fleld. Here e denotes the elementary charge and 

Q = 2g^5^|g;|, (16) 
/ 
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Fig. 1. Illustration of the Chiral Magnetic Effect in a very large homogeneous mag- 
netic field. The red arrows denote the direction of momentum, the blue arrows - the 
spin of the quarks. (1) Due to the very large magnetic field the up and down quarks 
are all in the lowest Landau level and can only move along the direction of the mag- 
netic field. Initially there are as many left-handed as right-handed quarks. (2) The 
quarks interact with a gauge configuration with nonzero Qw Assuming Q„ = — 1, 
this will convert a left-handed up/down quark into a right-handed up/down quark 
by reversing the direction of momentum. (3) The right-handed up quarks will move 
upwards, the right-handed down quarks will move downwards. A charge difference 
ofQ = 2e will be created between two sides of a plane perpendicular to the magnetic 
field. 

where g/ is the charge in units of e of a quark with flavor /. For Nf = 2 and 
Nf = 3 the relation above becomes Q = 2Qw and Q = ^Q^ respectively. An 
alternative derivation of this result is presented in Appendix B. 

Since Eq. ( !T6|) was obtained in the most ideal case, i.e. chiral limit and an 
extremely large magnetic field, Eq. (fT6l) has to be an upper limit. Therefore 
there is a maximum amount of charge that can be separated by a particular 
gauge field configuration in a homogeneous background magnetic field; in other 
words 



In our derivation we ignored the back-reaction due to the electric field 
created by the separating charges. We believe that this back-reaction can only 
give rise to a small suppression of the effect. This is because the size of the 
electric field will be much smaller than eB since there are very few particles 
involved. Moreover, in the physical case we are interested in, a heavy ion 
collision, the color forces will surely dominate. 

Since the separation is independent of color, the mechanism will not create 
a net color charge difference. Therefore we can safely ignore a gluonic back- 
reaction. 

4-2 Moderate Magnetic Field 

Now that we considered the ideal situation with an extremely large magnetic 
field, let us estimate the amount of charge separated in a moderate homoge- 




(17) 



/ 
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neous magnetic field. There are now three important scales in the problem, 
eB, the temperature T, and the size p of the configuration with nonzero Q^. 

In a moderate magnetic field not all spins will be aligned along B. Let us 
introduce the function (3f{uo) which denotes the degree of polarization along 
B of quarks with flavor / and energy smaller than uj. Since a magnetic field 
does not distinguish between left and right-handed particles, (3f{uj) will be 
the same for left and right-handed particles in thermodynamic equilibrium. 
Moreover, ignoring chemical potential, (3f{(jj) will be equal for particles and 
antiparticles. 

In a homogeneous magnetic field it is possible to calculate Pjiuj). The 
energy dispersion relation of a massless fermion is then given by 



u;=pi + 2\qfeB\n, (18) 

where n = 0, 1, 2, . . . denotes the Landau level. The n = level contains only 
one spin direction, all other levels contain both spin directions and are hence 
double degenerate. 

The number density of quarks in a magnetic field with energy smaller than 
u) reads 

N^ + N^ = r ^VzanOiu' - colMco,), (19) 

where a„ = 2 — 5nfi is the degeneracy factor and niup) = [exp{u!p/T) + 1]^^ 
is the Fermi-Dirac distribution function. The difference between the number 
density of fermions with spin parallel and antiparallel to the magnetic field 
yields 

I^T - Nl I = dpsOico' - pIMp,). (20) 

Hence the degree of polarization of quarks/antiquarks with energy smaller 
than uj is given by 

Let us try to find an approximation for Pf{uj). For 2\qfeB\ > u"^, Pf{uj) = 1 
since only the lowest Landau level contributes. For 2\qfeB\ = 0, Pf{uj) = 0. 
Therefore 

with 



/3/(^)^7(^^) (22) 



, X for X < 1 
l{x)={ - . (23) 

1 for X > 1 

We have computed Eq. ( l2Ti) numerically and found that Eq. ( l22i) is a reason- 
able approximation. 

At first hand, it might seem surprising that the temperature drops out 
the estimate for the degree of polarization. The reason for this is that we 
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Fig. 2. Crossectional view of a non-central heavy ion collision along the beam-axis 
(z-axis). The two nuclei have radii R, travel in opposite directions and collide with 
impact parameter b. The plane y = is called the reaction plane. The angle (p is 
an azimuthal angle with respect to the reaction plane. The region in which the two 
nuclei overlap contains the participants, the regions in which they do not contain 
the spectators. 

only considered the degree of polarization of a subset of all modes, namely 
the ones with energy smaller than u. If one would instead be interested in the 
polarization of the bulk of the particles, one should take ~ T, by which the 
temperature dependence naturally is recovered. 

Quarks with momenta much larger than the inverse size 1/p of the con- 
figuration with nonzero winding number will not be affected by these config- 
urations and hence do not change helicity. Therefore the quarks that change 
helicity will have momenta and hence energy smaller than ~ 1/p. Assuming 
that all quarks with energy smaller than 1/p are equally likely to change he- 
licity we find the following estimate for the expectation value of the amount 
of charge separated 



/ 

where $ = eBp^ is essentially the flux through a configuration with nonzero 
Qw For small fields {2\qfeB\ < p^) this reduces to 



5 Chiral Magnetic Effect in Heavy Ion Collisions 

Now that we know how much charge is separated by a single configuration 
with winding number Q^, let us estimate the amount of charge separated in 
a heavy ion collision. We will assume that hot deconfined chirally restored 
matter is created, so our results of the previous section will be applicable. 

As we have shown in the previous section, if eB ~ 1/p^ there is a good 
chance of getting charge separation. In a hot plasma, the typical size of 
a configuration with nonzero is of the order of the magnetic screening 
length, asT. The typical required magnetic fields are therefore of order a^T^, 




(24) 




(25) 



/ 
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which we estimate to be of order 10^ — 10^ MeV^, for the initial conditions 
typical at RHIC. In Appendix |X1 we show that such enormous magnetic fields 
can indeed be created in off-central heavy ion collisions. In the plasma, the 
magnetic field is to a good degree homogeneous and pointing in a direction 
perpendicular to the reaction plane = in Fig. ([2])) in the region of overlap 
of the colliding nuclei for the impact parameters of interest. 

We will now define and A''^ to be the total positive /negative charge 
in units of e above (a) and below (6) the reaction plane respectively. We have 
that iV^ = N^ + N^, where N"^ is the total positive /negative charge produced 
in a certain event. We define A± to be the difference between in charge on 
each side of the reaction plane A± = — N^. 

One generates an electric current along the magnetic field each time a 
transition from one vacuum to another is made. Let us now imagine such 
transition happening in the center of the hot matter. Locally a charge differ- 
ence will be created due to the current. But since the quarks creating this 
current will still have many interactions before they leave the hot matter, 
presumably transitions in the center due to the final-state randomization of 
quark momenta will not give rise to any observable charge separation. 

Now let us look to a transition near the surface of the hot matter. Sup- 
pose that the situation is such that the charge separation mechanism causes 
more positively charged quarks to move outwards and more negatively charged 
quarks the move inwards (the opposite situation is of course also possible with 
an equal probability). The positively charged quarks moving out of the sur- 
face will hadronize and can essentially freely propagate to the detector. On 
the other hand the negatively charged quarks will move inwards and still en- 
counter many interactions. As a result the azimuthal angles of the positively 
and negatively charged quarks will become uncorrelated. This is very similar 
to the observed jet quenching in heavy ion collisions. Only for very large im- 
pact parameters, when the size of the interacting region is about same as the 
scattering length, we expect correlations between positively and negatively 
charged particles to exist. 

Each time a transition is made both A+ and A_ will change. The expec- 
tation value of the change is either positive or negative with equal probability 
and given by 

±J2\qf\li2\qfm±M- (26) 

/ 

Here we used that the most probable transitions are those which connect 
two neighboring vacua, so Q„ = ±1. The transitions which simultaneously 
induce changes with Q„ > 1 are suppressed, therefore we have neglected 
them. We compute those transitions which are independent of one another, 
but taken over the time and volume of the system can induce multiple changes 
in the Chern-Simons charge. The functions C±{x±) are screening suppression 
functions in order to incorporate the effect that transitions near the surface 
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are much more likely to contribute to A±. We will take 



^±ix±) =exp{-\y±{x) -y\/X), 



(27) 



where A is the screening length. The functions y+{x) and y-{x) denote the 
upper and lower y coordinate of the overlap region. Hence \y±{x) — y\ is the 
distance from a point y to the upper/lower part of the surface. It follows from 
Fig. ([2D that 



y^[x) = -y-{x) = 



(x 



-6/2)2 -R + b/2<x<0 

^ L_ ^ " " , (28) 

v/i?2 - (x + 6/2)2 0<x <R-b/2 



here R denotes the radius of the nuclei and b is the impact parameter. Now 
one can easily infer from Eq. fl26l) that the screening works. It is equally 
likely to increase A± as to decrease it. However if for example y is near y^ 
the expectation value for increasing A± is much larger then the expectation 
value for decreasing it. Let us mention here that Eq. fl2^ is strictly speaking 
only valid for a constant homogeneous magnetic field. In the overlap region 
the magnetic field is to a good degree homogeneous around zero space-time 
rapditiy especially for large impact parameters. It has however a huge time 
dependence. So in reality there could be effects which are due to the rapid 
time variation of the field that we will ignore in this paper. 

Now we define to be the total number of raising/ lowering transi- 
tions. It of course holds that A^t = + ^t~- Furthermore we introduce 
At which is the difference between the number of raising and lowering tran- 
sitions, i.e. At = — . We can assume that all transitions happen 
independently from each other. Then the dynamics for At is exactly governed 
by a one- dimensional random walk. At each transition one has equal proba- 
bility to go up or down. Hence the expectation value of At will vanish. But 
since there are A^t transitions its variation is equal to v^iVt, so that 



dt 



d^x / dp 



dN, 



d-^x dt dp ' 



(29) 



where V denotes the volume in which the transitions take place. 

Now we can compute the variation of A^. Since we assume that all tran- 
sitions are independent from each other the total variation is the sum of the 
variation of all contributions. Hence it follows that 



(A|) 



n 2 



J2\QfhmQfeB\p' 
. f 



■ (30) 
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Likewise we can compute (A_|_A_). We find 



■ f 



■ (31) 



Let us remark here that if screening is absent we have (A^) = — (A_(_A_) since 
then ^± = 1. In the presence of screening however, ^+^- < |(^+ + so that 
KA+A_)|<(A|). 

Using the fact that p ~ {T"^ /as)~^^^ ~ l/{asT), and the expression for 
the sphaleron transition rate Eq. f|TT]) we can rewrite the last two formula's 
for small magnetic fields {2\qfeB\ < 1/p^) as 



d(A|) 
drj 



2k as 



-\ 2 



. / 



d(A+A_) 
d?7 



—4k as 



. / 



d^a;±^+(a;_L)^_(x_L) 



V, 



drr [e5(r,?7,x_L)]^ 
(32) 

7 

drr [e5(r,?7,x_L)]^ 
(33) 

where the proper time t = {t^ — z^y/'^ and the space-time rapidity rf = 
I log[(t + z)/{t — z)]. The volume integral is over the overlap region V± in the 
transverse plane. The time integral is from the initial time rj to a final time Tf. 
We have assumed here that the magnetic field does not change the sphaleron 
transition rate dramatically. We have also inserted a constant k which should 
be of order one, but with large uncertainties, since in our estimate of the scale 
size there is an uncertainty of order one which is taken to the fourth power. 
We assumed that the scale size factor in the sphaleron roughly canceled the 
factor of 385ci;| in the numerator, since the rate per unit volume should be 
of the order of the inverse size of the typical field configuration to the fourth 
power. 

At this point, we could take our exact expression for the magnetic field 
given in Appendix A and compute the rate. This we will do in a later paper for 
a variety of initial conditions. For now however we will use the approximate 
formulae for the magnetic field at the origin due to spectators and partici- 
pants. The approximations are given in Eqs. ( ]A.12[) and ( 1A.17I) . and valid 
for R/ sinh(yo) ^ ^ ^ -R- Adding the contributions of the spectators and 
participants we find 



eB ^ Za 



EM 



cexp(-Fo/2) 



i?l/2r3/2 



f{b/R)+Aexpi-2Yo)- 



(34) 



Here Z denotes the charge of the nucleus, R its radius and Yq the beam 
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rapidity. The function f{h/R) given in Eq. (lA.lSP and displayed in Fig. (1A.2P 
describes the impact parameter dependence of the participant contribution 
to the magnetic field to first approximation. The constant c = 0.59907... 
and aEM ~ 1/137 denotes the fine structure constant. Especially for large 
impact parameters the magnetic field at the origin is very homogeneous in 
the transverse plane within the overlap region, so that Eq. flMl) is a good 
approximation to the field at the surface of the interacting matter. For smaller 
impact parameters the magnetic field at the surface is somewhat smaller, so 
that we might overestimate the magnetic field at the surface using Eq. fl5^ 
for small impact parameters. 

Let us assume that the screening length A is constant in this short time. 
Then the integral over the transverse plane can be done since we used that 
the magnetic field is homogeneous to a certain degree. We define 

r r (35) 



h{b/R,X/R)^— I dVe+(^±)e-(^±). 



V, 



In the limit of A ^ oo these functions become equal to 



i5^^d^., = 2arcoos(4)-yi-j^. (36) 

We can now perform the integration over proper time. In the integration 
over proper time, the magnetic field is sharply falling. We can therefore to 
good approximation take = oo and find 



d(A|) 
drj 



. / 



2 

c^exp(-yo)/'(&/^)- 

Ti 



g{b/R,X/R) 



16 , , ,,.62^2 



+— cexp(-5/2yo)/(V^)^7^ + 4exp(-4Fc 



(37) 



We now must motivate a choice of initial time in the above equation. 
If the time at which the topological charge changing processes is less than 
the Lorentz contracted size of the system, then the magnetic field is at its 
maximum value, and our approximation of treating the spectator nucleons 
as being on an infinitesimally thin sheet breaks down. At the time when 
this approximation breaks down, the field should have achieved its maximum 
value. In this case, we should be able to approximate the initial time as 

Ti = (R e-^° (38) 

where C ~ 2. On the other hand, if the initial time for the topological charge 
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changing processes is larger than this, we should be able to approximate 



1/Q 



sat 



(39) 



where Qsat is the saturation momentum. The saturation momentum sets the 
time scale for the evolution of classical fields produced in the collisions. The 
choice between these two times scales is the maximum of the two possible 
values. 

For gold nuclei at 100 GeV per nucleon, both these choices are roughly 
the same and of the order of 0.1 — 0.2 fm/c. For copper at the same energy, 
the choice would be the inverse saturation momentum since a copper nucleus 
is much smaller than gold. For gold at lower energy, we would best choose 
the Lorentz contracted size scale, since the finite size of the nucleus becomes 
important. 

If we use the Lorentz contracted size scale, appropriate for gold, then 



d(A|) 
drj 



. / 



g{b/R,X/R) 



Similarly we find 



(40) 



d(A+A-) 



-AZ'^a^MK as 



. / 



h{b/R, \/R) 



c'nhiR)- 



(41) 



We expect that taking = 2 is reasonable. For two flavors we get (^^ g|)^ = 
25/81. Then we find for collisions at large impact parameters when screening 
is unimportant 



d{/\\) 100 



2 arccos 



h 

2R 



r' 



4i?2- 



(42) 



Now let us take b/R = 1.6 which roughly corresponds to 50 — 60% centrality. 
For gold-gold collisions {Z = 79 and E/A = 100 GeV ) we then find taking 
as = 1 and n = 1 

d{Al) 



drj 



0.1. 



(43) 



Let us stress again that this is a very rough estimate. For example taking the 
initial time twice as large will decrease our estimate by a factor 16. We could 
easily increase our estimate by an order of magnitude by taking k, = 10. 

In the next section we will connect (A^) and (A_|_A_) to some observables. 
Then we will discuss in detail the implications of Eqs. (HOll and (HTll . 



17 



6 Observables 



Due to the finite multiplicity of produced particles, one cannot observe charge 
separation in an individual event. The statistical fluctuations ~ will be 
much larger than the expected charge asymmetry. By taking an average over 
many events A± will of course vanish, since configurations with positive and 
negative Qw are very likely to be produced with equal probability. However 
can define correlators which are sensitive to the variation (A^). In this section 
we are going to connect these correlators to (A^). 

A set of very useful correlators for the study of charge separation was 
proposed by Voloshin [l9]. For each event one defines 



^ i=0 i=0 

where a, 6 = ± denotes the charge, N± is the total number of positively 
or negatively charged particles, and (pai denotes the azimuthal angle of an 
individual charged particle with respect to the reaction plane. 

Then in order to remove the multiplicity fluctuations one averages the 
correlators over A'e similar events. This averaging is called event mixing. In 
this way one can define the averaged correlators a , a++ and [49j, where 



flab = -^X] /(0a, 0b). (45) 

^ 71=1 

We feel that it is also useful to define correlators which are not divided by 

square of the total multiplicity of charged particles. Here we introduce h , 

and where 

bab = -^Y^^^^'^b^' (46) 



n=l 



with 



Na Nt 

The correlators can also be written as 

(7(0„0fe) = X,X,-y,n, (48) 

where 

X±^Y1 ('^n > ^± = E i"^^) ■ (49) 

i=l i=l 

An equivalent representation of the last two equations reads 

X±=/ d0cos(0)-f, r±= / d0sin(0)-f, (50) 
Jo d0 Jo d0 
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where dA^± denotes the number of charged particles in the azimuthal angle 
d0. Now one can see that Xj. essentially measures directed flow. By taking 
a symmetric interval in rapidity this quantity vanishes. On the other hand 
Y± measures the difference between charges on opposite sides of the reaction 
plane. Let us emphasize that if a measurement shows that a++/6++ is positive 
there really is on average more charge on one particular side of the reaction 
plane. Possible contamination by directed flow will always give a negative 
contribution to a_|_+/6++. The correlator a+_/6+_ measures correlations be- 
tween positively and negatively charged particles and the reaction plane. If 
such correlation exists and positively charged particles are emitted oppositely 
from negatively charged particles a+_/6+_ is necessarily negative. 

In order to relate the correlators to A± we have to make an assumption of 
the distribution of the charges in the azimuthal angle. If all charged particles 
affected by the Chiral Magnetic Effect would be really emitted perpendicular 
to the reaction plane in which case = 7r/2 or0 = 37r/2we would have 

Y± = A±. (51) 

As a result 

a++ = a__ = ^(A|), a+_ = ^^^(A+A_), (52) 

and 

6++ = 6__ = (A|), 6+_ = (A+A_). (53) 

where N± denotes the total number of charged particles in the corresponding 
T) interval. 

However, it might be more realistic that the particles will spread along the 
direction of the average magnetic field. If for example 

^"^'V. + iA.sinW, (54) 



we would find 



y± = ^A±. (55) 



1 7r2,.2> 1 TT^ 



In that case 
Similarly we find 



b^-,^b^^^-{Ai), 6+_ = -(A+A_). (57) 
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7 Discussion 



One should be aware that the observables defined in the previous section are 

V and T even. For vanishing net charge density per unit rapidity at = 0, 

the invariance under C predicts that a++ = a and 6++ = h Hence not 

only the Chiral Magnetic Effect contributes to the observables, but also some 

V even processes could do in principle. In order to find out whether the signal 
seen in Ref. [50] is really due to the Chiral Magnetic Effect one has to perform 
a detailed analysis and see whether calculations based on the Chiral Magnetic 
Effect are consistent with the data, and rule out more mundane explanations. 
Based on our main results Eqs. fHOj) and (14T|) we therefore will make some 
predictions. Let us make a cautionary statement here. We have made many 
approximations and there are many uncertainties. We expect that at best we 
might describe the data only semi-quantitatively. 

7.1 Absolute Value of Asymmetry 

In Eq. (1431) we estimated with large uncertainty (Aj_) = 0.1 at h/R = 1.6 
which corresponds roughly to 50-60 % centrality. Using that that there are 
about 24 charged particles produced within \7]\ < 0.5 (see Table [T]) we find 

a++ = a ~ 10""^. Now this estimate is very uncertain, we can easily increase 

or decrease it by an order of magnitude. However it seems comparing this 
number to the results presented in Ref. [50] that our prediction has about 
the same order of magnitude as seen in experiment. The main lesson from 
this is that although we are very uncertain about the value of the asymmetry, 
the Chiral Magnetic Effect does not give rise to an asymmetry which is many 
orders of magnitude away from the data obtained by experiment. 

7.2 Z Dependence 

The Z (charge/atomic number) is probably our most robust prediction. Clearly 
the magnetic field is proportional to Z. Since the polarization of the quarks is 
proportional to B, the correlators bab are proportional to Z^. If polarization 
would become maximal however, which we believe is unlikely, there will be 
no Z-dependence. In principle one could test Z-dependence by taking nuclei 
with the same A and different Z. This however could be difficult since it will 
require very accurate data. 

7.3 A Dependence 

The A dependence of our result depends upon which time scale is most im- 
portant, that associated with the Lorentz contracted size of the nuclei, or that 
of the inverse saturation momenta. The specific A dependence depends upon 
what range of impact parameter, beam energy and nuclear size we consider. 
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7.^ Beam Energy Dependence 

We are very uncertain about the beam energy dependence, since our result will 
depend strongly on what we take as our initial time. If the initial time scales 
with the longitudinal size of the nucleus, i?exp(— Fo), we expect no beam- 
energy dependence at all for bab- However if it turns out to be better to use the 
inverse saturation momentum for initial time, we can get a large dependence. 
In such case we always expect (as long as the pancake approximation stays 
good and a quark gluon plasma is formed) that and bat are smaller at 
larger beam energies, we never expect them to become larger. 

7.5 Centrality Dependence 

The centrality dependence of the correlators aab and bab has a strong depen- 
dence on initial time and could also have beam energy dependence. Let us 
assume that Xj = ^i?exp(— Fq), we can then use Eq. f HOj) to determine the 

centrality dependence of and b , b^^. For now we will take C = 2 

and A = 1 fm and present a plot of aab of our model for Gold-Gold. It could be 
very well argued that one should take other values, but these seem to us not so 
unreasonable. In order to compare our results with the data we will plot a++, 
a as a function of centrality and normalize the value in the 50-60 % central- 
ity bin to 7 X 10~^ which is the value found by experiment ^U\. To convert the 
impact parameter to centrality we used a table provided in Ref. {BT]- Using 
their results we also computed the number of charged particles produced in 
each centrality bin (see also Ref. [62] for experimental data). These results 

are summarized in Table [H Our result for a++ and a is displayed in Fig. [31 

Qualitatively the result agrees with the data presented in Ref. [5U] . 

7.6 Suppression of a^_ Relative to a^+ and a 

In Sec. [Sjwe found that |(A+A_)| < (A|) because of the screening effect. As 
a result |a+-| < a++. So it is very natural that the absolute value of a+_ is 
smaller than that of It easily follows from our results that 

|a+_| _ h{b/R,\/R) 
g{b/R,X/R)' 

where the functions g{b/ R, X/ R) and h{b/R,X/R) are given in Eq. fl35l) . We 
displayed |a+_|/a++ as a function of b/R in Fig. ([!]) for different values of 
X/R. To compute the functions g{b/ R, X/ R) and h{b/R,X/R) we used the 
screening suppression factor ^±{xj_) given in Eq. flTTI) . It could very well be 
that this simple suppression is to simplistic, and has to be modified somehow. 
Especially near the surface, and at large impact parameter, the function could 
be different. Nevertheless, one can infer from Fig. that for reasonable 
values of X/R one finds an order of magnitude suppression of a+_ compared 
to a++ for small impact parameters which is consistent with data presented in 
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(58) 



centrality 


h/R 


-^part 


(iV+ + iV_)/2 


- 5 % 


0.35 


344 


283 


5 - 10 % 


0.47 


317 


256 


10 - 20 % 


0.82 


222 


169 


20 - 30 % 


1.07 


153 


111 


30 - 40 % 


1.26 


102 


72 


40 - 50 % 


1.43 


64 


43 


50 - 60 % 


1.58 


38 


24 


60 - 70 % 


1.72 


20 


12 


70 - 80 % 


1.84 


9.5 


5.5 



Table 1 

Relation between centrality, impact parameter, number of participants A'part, and 
the average number of produced charged particles (A^+ + N^)/2 within |?7| < 0.5 for 
Gold-Gold collisions at 130 GeV per nucleon pair. 
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Fig. 3. Possible result of the Chiral Magnetic Effect for a++ and a as a function 

of centrality, for Gold-Gold collisions at 130 GeV per nucleon pair. 



Ref. ^SD]. For larger impact parameters the suppression is smaller, since the 
system size is in that case smaller. 

Since the screening length A is independent of the size of the system, we 
expect |a+_|/a_|_+ and to be larger for smaller systems. This implies 
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a++ 




Fig. 4. Possible result of the correlator |a-| \/a^ 

values of the screening length X/R. 



as a function of b/R for different 



that |a+_|/a++ and should be larger in Copper-Copper collisions 

than in Gold- Gold at the same centrality and energy. 



7. 7 Transverse Momentum Dependence 

We expect that the particles which make the helicity flip have momentum 
in the order of aT, hence we expect that the main contribution to our bab 
comes from particles which have low transverse momentum. However we can- 
not clearly determine what is low. To make a safe prediction we expect the 
main contribution to bab come from particles which have transverse momentum 
smaller than 1 GeV. 



7.8 Particle Species Dependence 

We think the main contribution to bab should come from pions, since the Chi- 
ral Magnetic Effect also creates an isospin difference, but no baryon number 
difference if Nf = 2 (see also Appendix B). However at this point we cannot 
make definite statements about the ratios protons to pions or kaons to pions 
that contributes to bab- 

7.9 Order Parameter for Chiral Symmetry/Confinement? 

Since it is absolutely necessary for the Chiral Magnetic Effect to operate that 
chiral symmetry is restored and matter is in the deconfined phase, measure- 
ment of aab/bab if due to the chiral magnetic effect could also serve as an order 
parameter for chiral symmetry and confinement. If the temperature (beam 
energy) is too low such that there is chirally broken confined matter created 
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we think it is likely that aab and bat should vanish or at least become much 
smaller. 

8 Summary 

The Chiral Magnetic Effect allows one to in principle probe the effects of 
topological charge changing processes in QCD. These phenomena are of great 
interest in field theory, and of great potential application in cosmology. In 
QCD, such phenomena are perhaps at the generation of nucleon mass, and 
also are associated with dynamical breaking of V and CV symmetries. 

A confirmed observation of such an effect would be of great import. A 
confirmed observation must also be subject to alternative hypothesis. Our ex- 
planation has the advantage that many of its features may be varied by choice 
of beam energy, and the Z and A of the colliding nuclei. The computations 
we present are at present semi-quantitative estimates, that can be refined by 
more careful treatment of the magnetic field, and a better computation of the 
underlying topological charge changing processes. 
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A Magnetic Field in Heavy Ion Collisions 

Let us estimate the size of the magnetic field achieved in a heavy ion collision. 
We consider two similar nuclei with charge Z and radius R. These nuclei are 
traveling with rapidity Yq in the positive and negative ^-direction. At t = 
they collide with impact parameter b. We take the center of the nuclei at t = 
at X = ±6/2 so that the impact parameter direction lies along the We 
depicted this situation in Fig. [2l 

The nuclei in typical heavy-ion collision experiments are nearly travel- 
ing with the speed of light. For example gold-gold collisions at RHIC are 
performed at center of mass energies of 200 GeV per nucleon pair which cor- 
responds to Yq ^ 5.4. In that case the Lorentz contraction factor 7 is equal to 
100. Hence the nuclei are Lorentz contracted in the z-direction to 1 percent of 
their original size. Therefore to good approximation we can consider the two 
including nuclei as pancake shaped. As a result we can assume that number 
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densities of the two nuclei are given by 



''^^^'^^ ^ |^V^^^^^<±V2j^^±(a.'J (A.l) 

here x'^ denotes a vector transverse to the beam axis and 

^±(cr'J=^[it:2-(a;'^±6/2f], (A.2) 

The functions 6'±(cc^) are the projections of the nuclei on the plane transverse 
to the beam axis and will also be used to separate the spectators from the 
participants. The number densities are normalized such that 

y A^x'^p^{x'^) = 1. (A.3) 

First let us quote the magnetic field at position x = {x±,z) caused by 
a particle with charge Z moving in the positive z-direction with rapidity Y. 
At t — the particle can be found at position x'j_. This magnetic field can 
be obtained by either boosting the electric field of a charge Z or using the 
Lienard-Wiechert potentials. One finds {e — |e|) 

eB{x) = ZaEMsmh{Y) (a;^^ - x^) x 

[{x'^-x^y + itsinhY-zcoshVyf^ 

Now we can estimate the strength of the magnetic field during a heavy ion 
collision. We are only interested in t > 0, i.e. just after the collision. Then we 
can split the magnetic field in the following way 

B = Bt + B: + B; + B;, (A.5) 

where Bf and B^ are the contributions of the spectators and the participants 
moving in the positive or negative ^-direction respectively. We assume that 
the spectators do not scatter at all and stay traveling with the beam rapidity 
Yq. It follows that {Yq > 0) 

eB^{T,r),x^) = ±ZaEMsmh{YoTr]) J d^x'^p±{x'^)[l - ^^«)] 

[(x'^-Xj_)^ + T^smh{YoTv)^f^' 

where the proper time t — {t"^ — z^Y^^ and the space-time rapidity 77 = 
'^log[{t + z)/{t-z)]. 

Many particles will be created by the interactions between the participants. 
However, the number produced positively and negatively charged particles will 
be nearly equal. Moreover the expansion of these produced particles is almost 
spherical. We therefore expect that the contribution of the produced particles 
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to the magnetic field is very small and will neglect this contribution. Hence 
we only have to take into account the contribution of the participants which 
were originally there. It is know that these participants stay traveling along 
the beam axis according to the following normalized distribution 



fiY) 



2 sinh(aFo) 



-Fn < y < Fn. 



(A.7) 



Experimental data shows that a ~ 1/2, consistent with the baryon junction 
stopping mechanism (see p3] and references therein). 

The contribution of the participants to the magnetic field is hence given 

by 



^0 



eB^{T,r],x^) = ±ZaEM / d^^l / dY f{Y) smh{Y t v)p±{O0Ti^'±)x 



Yo 



[{x'j_ — x_i_y + sinh(F =f 77) 



2l3/2 



. (A.8) 



We have evaluated the magnetic field numerically at the origin (r/ = 
0, a;_L = 0) in which case it is pointing in the y-direction. We took collid- 
ing gold ions [Z = 79, R = 7 fm) with different beam rapidities {Y = 4.19 
and Y = 5.36). The results are displayed in Figs. lA.ll and IA.2I Clearly 
enormous magnetic fields are created in non-central heavy ion collisions. 
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Fig. A.l. Magnetic field at the center of a gold-gold collision, for different impact 
parameters. Here the center of mass energy is 62 GeV per nucleon pair (Yq = 4.2). 

For our purposes it is useful to have some analytic expression of the mag- 
netic field. We can find reasonable approximations for r > 2R/ sinh(yo)- First 
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Fig. A. 2. Magnetic field at the center of a gold-gold collision, for different impact 
parameters. Here the center of mass energy is 200 GeV per nucleon pair (Yq = 5.4). 

we will consider the spectators, then we will discuss an approximation for the 
participants. We will perform both approximations at the origin {x± = 
and 1] = 0). In that case the magnetic field is pointing in the y-direction, 
eB = eBSy. Especially for large impact parameters the magnetic field at 
the origin will be a good estimate for the magnetic field at the surface of the 
interacting region, since the magnetic field in the overlap region is to a good 
degree homogeneous in the transverse plane. 

A.l Spectator Contribution for r > R/ sinh(yo) 

For r > i?sinh(yo) the denominator of the integrand of the spectator contri- 
bution Eq. flA.6p can be approximated by sinh(yo)'^- Hence we find 



where 



with 



AR 

eBs ^ ZaEM (iy^'9{-'2YQ)—g{h/ R), 



g{h/R) = Y,9±{h/R), 

± 



(A.9) 
(A.IO) 



g±{h/R) = T^I dV^p±(a.l)(l - e^{x'^))x'. (A.ll) 
We find that to very good approximation g±{b/R) = b/R. As a result 
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eBs ^ ZaEMexp{-2Yo) — . 



(A.12) 



27 



A. 2 Participant Contribution for R/ sinh(yo) ^ ^ R 

At midrapidity {rj = 0) the participant contribution to the magnetic field 
Eq. (lA.SIl contains the following integral over the beam rapidities 



2 sinh(aFo) J^Yo [/?^ + smh^Y)]^^ 

where = {x'^ — x±Y /t"^. We are interested in finding an approximation to 
Eq. ( 1A.13[) for 1 < /3 < sinh(yo)- The first step is to write 



exp(aFo) 2r3 [/?2 + i exp(2F)]3/2 ' 

The last integral can be computed exactly. However, since it is quite compli- 
cated for a = 1/2 and because we are anyway interested in the limit where 
1 < 13 < sinh(Yo) we will expand the result in /?/ sinh(yo)- This gives for 
a < 1, 

I ^ ce-xY>{-aYQ)-—— (A.15) 

where 

2^+"ar(l-a/2)r(3/2 + a/2) 
(1 + ajVTT 
For a = 1 /2 we find c = 0.59907 . . .. 

Now let us try to approximate the magnetic field at the origin (77 = 0, x^ = 
0) using the approximation for /. Taking a = 1/2 we find that 

eBp ^ cZaEMexp{-Yo/2)^^^j^fib/R), (A.17) 

where 

f{b/R) = J2f±ib/R), (A.18) 

± 

with 

f^{b/R) = tR'" I d'x'^p^ix'MOj^,- (A.19) 
We have computed f{b/R) numerically and displayed the result in Fig. flA.2p . 



B Alternative derivation of Chiral Magnetic Effect 

We will show in a somewhat more formal way than in Sec. |4]that a homoge- 
neous background magnetic field in the presence of a gauge field configuration 
with nonzero Qw gives rise to a vector current = [if)'^ ^ip) a along the direc- 
tion of the magnetic field. 

We take the background magnetic field to point in the z-direction so that 
only the following component is non-vanishing, 

is = (^73^')^ = (^icrs^R)^ - (V'Ic^3^l)a- (B.l) 
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Fig. A. 3. Plot of the function f{b/R). This function describes the impact parameter 
dependence of the approximated participant contribution to the magnetic field and 
is defined in Eq. (|Al^ . 

Initially, there is no difference between the number of left-handed and right- 
handed handed modes, therefore js vanishes and there is no charge separation. 
Clearly Eq. ( IB.ll) shows that if left- and right-handed modes are polarized in 
the same direction (as it should in a large magnetic field), a charged current 
will only form if there is a difference between and A^^. More explicitly in 
the limit of a very large magnetic field {eB ^ or eB ^ where Pf is 
the Fermi momentum) 

J d^a; {^]i^L(^3^R,L)A = sgn(g) J d^x {iP^r^l'^^ r^l) a- (B.2) 

Here q denotes the charge of the particle and reflects the fact that particles 
with different charge are polarized oppositely in a magnetic field. We now 
immediately find that in the limit of large magnetic field the total current for 
a flavor / is given by 

4 = j d'xjiix) = sgn(g;)(iV^ - iV/). (B.3) 

In the chiral limit this implies that at t = oo an electromagnetic (EM) and 
a baryon (B) current is formed which magnitude equals 

J3^^ = 2g„ I^/I' = E sgn(g^). (B.4) 

/ / 

If we assume that all quarks were located within a finite volume, this 
electromagnetic current generates a charge difference Q between two sides of 
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a plane perpendicular to the magnetic field of size 

/ 

This difference is the same as we found before in Eq. (fT6|) . 

In the same way we find the following baryon number difference 

i? = 2^j/ = |g^5]sgn(g;). (B.6) 
/ / 

The baryon number difference vanishes for Nf = 2, while for Nf = 3, B = 
IQw Since B times the average charge of a baryon is smaller than Q, we 
expect the charge difference mainly to be realized in mesons, like tt"*" and vr". 

In this appendix we considered a gauge field configuration with nonzero 
Q„ in the presence of a very large homogeneous magnetic field. It would 
be interesting to extend these results by studying the behavior of the vector 
current in a smaller magnetic field. 
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